Surface criticality in random field magnets 
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The boundary-induced scaling of three-dimensional random field Ising magnets is investigated close 
to the bulk critical point by exact combinatorial optimization methods. We measure several expo- 
nents describing surface criticality: /3i for the surface layer magnetization and the surface excess 
exponents for the magnetization and the specific heat, /3s and Us. The latter ones are related to the 
bulk phase transition by the same scaling laws as in pure systems, but only with the same violation 
of hyperscaling exponent 6 as in the bulk. The boundary disorders faster than the bulk, and the 
experimental and theoretical implications are discussed. 
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I. INTRODUCTION 

The presence of quenched randomness leads to many 
differences in the statistical behavior if compared to 
"pure systems" . This is true in many phenomena as 
transport properties in, for instance, superconductors, or 
in a rather wide range of cases in magnetism. Consider 
a domain wall in a magnet, which gets pinned due to 
impurities. The scenario may vary according to the sym- 
metries of the system and to the character of the disorder, 
but is described, in most general terms, by an "energy 
landscape" which develops a rich structure due to the the 
presence of pinning defects [l| . 

The most usual and convenient example of such mag- 
nets is given by the Ising model -universality class. Disor- 
der is normally introduced as frozen random bond" and 
"random field" impurities, which can change dramati- 
cally the nature of the phases of the model and the char- 
acter of the phase transition. Strong enough bond dis- 
order creates a spin glass -state, while the random fields 
couple directly to the order parameter, the magnetiza- 
tion. 

The criticality in such models is usually studied by 
finite size scaling, to extract the thermodynamic behav- 
ior. However, real (experimental) systems are finite and 
have boundaries. These break the translational invari- 
ance and create differences in the critical behavior be- 
tween the boundary region and the bulk. The related 
phenomenon is called "surface criticality", and essential 
is that a whole set of new critical exponents arises, to 
describe the behavior of various quantities at and close 
to surfaces 0, |^. Here, we investigate by scaling ar- 
guments and exact numerical methods this phenomenon 
in the case of the random field Ising model (RFIM), in 
three dimensions (3d). In this case, the RFIM has a bulk 
phase transition separating ferromagnetic and paramag- 
netic states. 

The central question that we want to tackle is: how 
do disorder and the presence of boundaries combine, in 
a system where the critical bulk properties are already 
different from pure systems? Though disordered mag- 



nets have been investigated earlier for the case of weak 
bond-disorder i^,, ^] , both spin-glasses - a possible future 
extension of our work - and the RFIM have not been 
studied j^. One general problem of the 3d RFIM has 
been how to observe the critical behavior, and under- 
standing the boundary critical behavior provides an in- 
dependent, novel avenue for such purposes Such 
experiments are done on a number of systems from di- 
luted antifcrromagnets in a field, [tI B| ) to binary liquids 
in porous media, and to relaxor ferroelectrics 



The particular characteristics of the RFIM is a com- 
plicated energy landscape, which manifests itself e.g. in 
the violation of the usual hyperscaling relation of ther- 
modynamics, and in the existence of an associated viola- 
tion exponent 6 and several consequences thereof. This 
is analogous to, for instance, spin glasses, and further- 
more for surface criticality presents the question how the 
broken translational invariance combines with the energy 
scaling. Our results imply that this can be understood 
by scalings that include both the bulk correlation length 
exponent and the bulk 6 and novel surface exponents. 
Moreover, though the bulk RFIM 3d phase transition 
has been notoriously difficult experimentally, the bound- 
ary order parameter, say, should be quite sensitive to the 
control one (temperature, in experiments and disorder 
here) and promises thus to make the surface criticality 
experimentally observable. 



In the next section we overview the theoretical picture, 
as applied to the RFIM. Section 3 presents the numerical 
results, where the emphasis is two-fold. We discuss the 
surface criticality on one hand, and on the other hand 
the decay of a surface field induced perturbation is ana- 
lyzed, since it has characteristics peculiar to a disordered 
magnet, in contrast to pure systems. Finally, Section 4 
finishes the paper with a discussion of the results and 
future prospects. 
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II. SURFACE CRITICALITY 

The RFIM Hamiltonian with a free surface 5" reads 
HnFiM = -J ^ (JiUj-Ji ^ (TiO-j-^/iiO-i, (1) 

where J is the bulk (nearest neighbour) interaction 
strength while Ji describes the strength of the surface in- 
teraction, in general different from J. ai take the values 
±1. For simplicity, the random fields hi obey a Gaussian 



probability distribution P{hi) = 



/2ttA 



exp 



One might 



with a zero mean and standard deviation A. 
have also external fields such as a bulk magnetic field h 
and a surface magnetic field hi at S. 

Being governed by a zero temperature fixed point, the 
phase transition of the 3d RFIM can also be studied at 
T = 0, where it takes place at a critical A^. The transi- 
tion is of second order though it also exhibits some first- 
order characteristics: the order parameter exponent f3 is 
very close to zero Elllllli- The surface criticality of 
the 3d RFIM is simplified by the fact that the lower crit- 
ical dimension is two 0, ^| , thus in the absence of a 
surface magnetic field hi just an ordinary transition can 
take place. The surface orders only because the bulk does 
so, and the transition point is the bulk critical point. 

Even in this case, there is a wide variaty of surface 
quantities. Derivatives of the surface free energy fs (sur- 
face ground state energy at T = 0) with respect to 
surface fields, as the surface magnetic field hi, yield 
local quantities (e.g. the surface layer magnetization 
nil = —dfs/dhi), while derivatives of fs with respect to 
bulk fields produce excess quantities, such as the excess 
magnetization nis — —dfs/dh, defined by 



If S 
— / d'^x r7i(x) = mb+ y^T-s + 0{L~ 



(2) 



where m(x) is the (coarse grained) magnetization at x 
and V r^L"^ and S are the sample volume and its surface 
area, respectively. One also obtains mixed quantities by 
taking second or higher derivatives of fs. We focus on the 
critical behavior of the local and the excess magnetization 
(mi and m^) as well as the excess specific heat Cs- 

The RFIM bulk critical exponents are related via the 
usual thermodynamic scaling relations, see Table The 
hyperscaling relations, however, have the modified form 



2-a^v{d-e), 



(3) 



with the additional exponent 9 16]. The usual way to 
relate the surface excess exponents to bulk exponents is 
to note that from the conventional hyperscaling (Eq. Q 
with 9 = 0) it follows that the singular part of the bulk 

free energy /^^^'"^•' scales with the correlation length ^ as 
^ gy making the analogous assumption for 



the surface free energy, /; 

Us = a + 



(sing) 



one finds 0] 
(4) 



In the case of the RFIM the above becomes less clear: 
does the ^-exponent get modifed? We assume that the 
exponent 9' in j]''™^' ~ ^-{d-i-e ) j^^^y -^^ general be 
different from the bulk exponent 9, and obtain 

as = a + i^-i^{9 - 9'), 
(is = P-v + v{9-9'). 

To derive Eq. 0, the scaling form 



(5) 
(6) 



t^ "'^Eslh/Jt (t+Z')] is used for the singular part of the 
excess ground state energy density E'f™^^ (which takes 
the role of the excess free energy at T = 0), with 
t = (A — Ac)/ J, Eq. (O and the Rushbrooke scahng 
law a + 2P + "f = 2. 7 is the exponent describing the 
critical behavior of the bulk susceptibility. Scaling rela- 
tions relating /3i to other 'local' surface exponents can 
also be derived, but it cannot be expressed in terms of 
bulk exponents alone. 



Quantity 


Definition 


Exponent 


excess magnetization 


nis = 


dfs 

dh 


TUs ^ 


-(- 


-tr^ 


excess specific heat 


Cs = 


d^fs 
dAdJ 


Cs - 




— as 


surface magnetization 


mi = 


dfs 
dh-i 


mi 


-(- 





TABLE I: Surface quantities in terms of the surface free 
energy fs, and the corresponding critical exponents (t = 
(A — Ac)/ J). Note that T = so that one uses instead 
of a free energy the ground state energy. 



III. NUMERICAL RESULTS 

The exact ground state (GS) calculations are based on 
the equivalence of the T — RFIM with the maximum 
fiow problem in a graph fl7l |: we use a polynomial push- 
relabel prefiow-type algorithm [1^.19,] . If not stated oth- 
erwise, we study cubic systems of size L^, L < 100. Free 
boundary conditions are used in one direction (the free 
surface under study) while in the remaining ones periodic 
boundary conditions are imposed. The maximal statisti- 
cal error in what follows is of the order of the symbol size 
used, so the error bars are omitted. Note that since in the 
present case only the ordinary transition is possible, the 
critical exponents should be independent of the surface 
interaction Ji. Complications arise, however, since in 2d 
the RFIM is effectively ferromagnetic below the break-up 
length scale Lh, which scales as Li, ~ exp [A{J/A)^] (see 
Fig. ^1 [20I I2H . This means that the surfaces have a 
tendency to be ordered "an sich" , and to see the true or- 
dinary transition behavior, one needs L > L^. Thus, we 
use substantially weakened surface interactions Ji <C J 
to circumvent this problem. 
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FIG. 1: The break-up length scale Lt of the 2d surface layer of 
the 3d RFIM with a strongly paramagnetic bulk, J — 0.05A, 
vs (Ji/A)^. Lb is estimated by looking for a value of Ji such 
that the surface will be totally ordered with probability 1/2 
while keeping A and L fixed. The solid line corresponds to 
A = 2.1. 



Using the bulk value v = 1.37 ± 0.09 T^l, one obtains 

/3i = 0.23 ± 0.03. (9) 

Fig. ^depicts miL^^^" versus (A - Ac)i^/'', and with 
Pi/iy = 0.17, ly = 1.37 and A^/J = 2.27 one indeed 
obtains a decent data collapse. With Ji « J, however, 
plotting mi (Ac) versus L produces a slightly different 
exponent, Pi/v « 0.15, and we could not get good data 
collapses, probably due to the fact that Li, is large. 
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FIG. 2: Mean absolute value of the surface layer magnetiza- 
tion mi as a function of A/ J for various L, Ji = J. The 
dashed vertical line corresponds to the critical point of the 
infinite system, A/J — 2.27. 



A. Surface layer magnetization 

Fig. 121 shows an example of the magnetization mi of 
the surface layer close to Ac, obtained directly from the 
spin structure of the GS. We assume the finite size scaling 
ansatz 



(7) 



where mi is a scaling function. At the critical point A — 
Ac, Eq. O reduces to mi L^^^^". Fig. Olis a double 
logarithmic plot of mi versus L at Ac/ J — 2.27 for three 
Ji-values. All three are consistent with 



^i/i^ = 0.17 ±0.01. 



(8) 
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FIG. 3: A log-log plot of the surface layer magnetization mi 
as a function of the system size L at criticality, A/J — 2.27, 
for various Ji/J <C 1. The solid lines depict fits, with f5\/v = 
0.17 ± 0.01 for all three cases shown 
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FIG. 4: A scaling plot of the surface layer magnetization mi 
in the case Ji = 0, J = 1, using Ac = 2.27, u = 1.37 and 
13-i = 0.23. 



B. Surface excess magnetization 

For the surface excess magnetization TOs, we use the 
finite size scaling ansatz 



^L-P^'''m,[{A~Ac)L^'''l 



(10) 
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where rhs is a scaling function. Since Pi was found to be 
independent of Ji/J as long as Ji/J ^ 1 (in the limit 
L — !■ cxD, the independence of the exponents on Ji/J 
should hold for any Ji/J), one expects the same to ap- 
ply for the other exponents as well and we thus consider 
here only the case Ji/J = 0.1. At the critical point, nis 
grows almost linearly with L (Fig. EJ, with the expo- 
nent —f3s/v = 0.99 ± 0.02. This yields, by again using 
V = 1.37 ±0.09, 

/3, = -1.4 ±0.1. (11) 
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FIG. 5: A log-log plot of the excess magnetization as a 
function of the system size L for A/ J = 2.27, Ji/J = 0.1. 
A background term of magnitude 1.07 has been substracted 
from TTis to see the power-law behavior. The solid line is a 
power- law fit, with —Ps/f = 0.99. 



C. Surface specific heat 

In GS calculations, the specific heat is computed (recall 
T = 0) by replacing the second derivative of the free 
energy / with respect to the temperature by the second 
derivative of the GS energy density E with respect to 
A or J dE/dJ is the the bond part of E, Ej = 

L^'^J2{ij) <^i<^j- The excess specific heat exponent as is 
estimated according to Ref. (where the bulk one was 
considered) . The singular part of the excess specific heat 
obeys 

C^/'"s) = L"=/'^C'J(A - Ac)Li/''], (12) 

from which by integration it follows for the singular part 
of the excess bond energy at criticality, 

47^' (L, A = Ae) = ci + c2L("-l)/^ (13) 

where ci and C2 are constants. Fig. is a plot of the 
excess bond energy, with Ji/J = 0.1, at the bulk critical 



point. The fit using Eq. (|13|l results in {ag — ^)/v = 
0.22 ± 0.03, corresponding to 

= 1.30 ±0.05. (14) 
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FIG. 6: A plot of the absolute value of the excess bond energy 
Ej,s as a function of L for A/J = 2.27, Ji/J = 0.1. The 
solid line corresponds to a fit of the form of Eq. I|18|l . with 
ci = 1.1292, C2 = 0.9756 and [as -l)/v = 0.22. 



D. Magnetization decay close to the surface 

Finally we discuss the behavior of the magnetization 
profiles m{z) (i.e. magnetization as a function of the dis- 
tance z from the surface) , in the case the spin orientation 
at the surface layer is fixed. This corresponds to applying 
a strong surface field hi . These are of interest as they re- 
flect spin-spin correlations close to the surface, as studied 
in Ref. [23| in the slightly different context of comparing 
two replicas with opposite hi. For the RFIM close to 
the infinite system bulk critical point, 'm{z) is affected 
by the fact that for numerically feasible system sizes the 
bulk magnetization is close to unity and decreases very 
slowly with increasing system size (due to the small value 
of 13) This is demonstrated in the inset of Fig. [7| 

where the distribution of bulk magnetization m;, at the 
critical point can be seen to be strongly peaked around 
mil = ±1. 

One can now distinguish three scenarios from sample 
to sample: if |mb| « 1 the applied strong surface field hi 
may have the same or opposite orientation, or finally the 
bulk magnetization mf, may be close to zero. In the first 
case, the hi induced spin configuration will be close to 
the one in the absence of the field. In the second case, hi 
will either force rrib to change sign altogether (producing 
again a flat profile with m(z) « ±1) or induce an inter- 
face between the two regions of opposite magnetization, 
as in Fig. [71 The third one has a small probability, and 
thus will not contribute much to the ensemble averaged 
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magnetization profile. The average magnetization pro- 
file {m{z)) can then (for a finite system, at the infinite 
system critical point) be well approximated by writing 

{m{z)) ^a + b{m,f{z)). (15) 

Here a and b are weight factors, here constant but in 
general function(s) of L, that tell the relative weight of 
samples where the magnetization changes inside due to 
the hi. 

{mif{z)) ^ J dwdzoPwiw)Pzoizo)m{z,zo,w) (16) 

is the profile one would obtain by averaging only over 
"single sample" profiles to(z, zq, w), corresponding to an 
interface of width w and position zq (with probability dis- 
tributions Pw and Pzo , respectively) . A simplified model 
for m{z, zo, w) is shown in Fig. |S| 

From the exact ground state calculations, we identify 
the profiles corresponding to such interface configura- 
tions. This is done by demanding that such profiles have 
a region where m{z) < —0.9 (when hi ^ 0). The inter- 
face width is defined as w = Z2 — zi, where zi and Z2 are 
the smallest z's such that m{zi) < 0.9 and m{z2) < —0.9, 
respectively. The interface position zq is then given by 
Zq = {zi + Z2)/2. By counting the fraction of such pro- 
files, we can estimate a and b in Eq. (|15|l . These have 
the approximate values of 0.39 and 0.61, respectively (for 
a system of size 40x40x80). By using Eqs. ((TK)l and ltTB|) 
with m(z,zo,w) presented in Fig. |H1 as well as the dis- 
tributions Pw and Pz„ measured from the ground state 
calculations, one indeed obtains an average profile {m(z)) 
that is in reasonable agreement with the true one, see Fig. 

El 

The average magnetization profile {m(z)) decays 
slowly with the distance z, not quite reaching zero at 
the opposite edge of the system in the case at hand. 
However, a typical value of m{z) will be close to ±1 
for all z, which persists for accessible system sizes due 
again to the small value of /3. One may thus observe 
effects reminiscent of violation of self-averaging, and 
this would be true also if one would measure the av- 
eraged difference (|m(z) — mcs{z)\) between the field- 
perturbed and GS configurations, and the higher mo- 
ments thereof. These results illustrate simply how the 
quasi-ferromagnetic character of the 3d RFIM ground- 
state influences such perturbation studies, a consequence 
of the in practice limited system sizes one can access in 
simulations. 



IV. CONCLUSIONS 

In this work we have studied with combinatorial op- 
timization and scaling arguments surface criticality in a 
random magnet, the 3d RFIM. The surface layer magne- 
tization exponent (3i is more than an order of magnitude 




FIG. 7: Main figure: A typical example of a magnetization 
profile, taken from a single sample, where due to a strong 
positive surface field /ii at 2; = an interface has formed 
between two regions of opposite magnetization. Inset: Distri- 
bution of the bulk magnetization mt with periodic boundary 
conditions, 2000 samples. A/J = 2.27, system size 40x40x80. 
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FIG. 8: A simple model for a single-sample magnetization 
profile m{z,ZQ,w). The interface is characterized by the pa- 
rameters position Zq and width w. 



larger than the extremely small bulk value [T^ 113, Hsf . 
Experimentalists have reported much larger values for /3 
0, H, 1 which in fact are rather close to our estimate 
for Pi . An intriguing possibility in this respect is the di- 
rect observation of the surface order parameter in relaxor 
ferroelectrics via piezoelectric force microscopy p^ . 

The excess exponents as and (3s, when inserted into 
the scaling relations ||SJ) and both yield very small 
values for the correction term 1^(9 — 9'), assuming a « 0, 
(3 w 0.02 and i/ « 1.37 1^. This suggests that in fact 
9' = 9, and the excess exponents are related to bulk 
exponents by the usual scaling laws valid for pure sys- 
tems, Eq. 101). The numerically obtained description 
of the ordinary surface transition uses the bulk correla- 
tion length exponent as in pure systems. All this would 
merit further theoretical considerations and could also be 
checked in the four-dimensional RFIM (2^ , whose phase 
diagram is also more complex due to the 3d surfaces 
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FIG. 9: Main figure: A comparison between the numerical 
{m{z)) (sofid line, averaged over 3000 samples) and that ob- 
tained by using Eqs. 11511 and HIGH with m{z, zo,w) as in Fig. 
|H|(dashed line). Inset: Distributions of the interface position 
Pzo{zo) (solid line) and width Pu,{w) (dashed line) obtained 
from the simulations. A/J — 2.27, system size 40x40x80. 



which have independently phase transitions. The spin- 
spin correlations close to the surface and the magnetiza- 
tion profiles in the presence of boundary perturbations 
have been studied, similarly to the context of looking for 



self-averaging violations [2J|. It would be interesting to 
investigate this aspect in more detail, but in our numerics 
the most transparent features are due to the two-peaked 
magnetization distribution of the groundstates, without 
a perturbing field. 

On a final note, the observations here concerning sur- 
face criticality in a disordered magnet - with a compli- 
cated energy landscape - extend directly for instance to 
spin glasses to a wide class of non-equilibrium 

systems (see 27;], also for experimental suggestions). 
Two evident possibilities are looking for the same phe- 
nomenology in 3d Ising spin glasses, and in the 3d zero- 
temperature non-equilibrium RFIM. In the former case, 
the free surface of a system at T > is in analogy to 
the zero temperature 3d RFIM case inherently disordered 
(the 2d spin glass has a, T ~ phase transition). In the 
second case, the situation is much more akin to the one 
at hand ([23) and one should consider as the order pa- 
rameter the remanent surface magnetization after a de- 
magnetization procedure. 
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